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Source  parameters  of  earthquake  of  31  ? 3.5  have  been 
analyzed  usiny  data  fro  ’ the  analog  and  digital  outputs  of  the 
HG-LP  stations.  Mosul to  fre  these  stations  are  co  pared  v.'itk 
those  of  the  7,7 o 07  stations. 


INTRODUCTION 


Since  1970  eleven  stations  with  high-gain  (magnifications 
up  to  150.000)  long-period  (peak  between  30  and  50  s)  systems  re- 
cording in  analog  and  digital  form  have  been  put  in  operation  (Sa- 
vino  et  al. , 1972).  These  stations,  that  will  be  refered  to  as 
HGLP,  have  been  used  in  the  determination  of  propagation  and  source 
mechanism  parameters,  comparing  their  performance  with  that  of  the 
VA7SSN  stations.  The  parameters  selected  in  this  study  ore  for  the 
propagation  characteristics  group  and  phase  velocities  of  surface 
waves  and  attenuation  and  for  the  source,  the  dimensions  of  the 
fracture  plane,  seismic  moment  and  the  stress  conditions.  Because 
of  the  response  of  the  HGLP  systems  which  peaks  at  about  30  sec, 
special  emphasis  has  been  placed  on  the  determination  of  the  about 
mentioned  parameters  by  means  of  surface  wove s. 

A very  important  step  in  the  analysis  o^  the  HGLP  data 
is  that  of  the  preliminary  processing  which  involves  analog  to 
digital  conversion,  filtering  and  spectral  analysis.  Reduction  to 
ground  motion  requires  also  knowledge  of  the  response  characteris- 
tics of  the  systems.  The  HGLP  instruments  hove  two  outputs,  or.e 
analog  (HGLP-AUL)  with  conventional  photographic  recording  the 
other  digital  (HGLP-DIG)  on  magnetic  tape.  Digitalization  of  the 
analog  records  have  been  done  from  enlarged  copies  of  the  microfi 1 - 


In  this  respect  the  data  from  earthquakes  relatively  large  and  near 
to  the  station  are  practically  impossible  to  analyze.  Amplitudes 
are  too  large,  lines  fainted,  and  traces  entangled.  The  digital 
outputs  should  solve  this  problem,  but  it  creates  some  of  its  own. 

In  order  to  study  dhe  capabilities  of  the  digital  outputs 
of  the  HGLP  stations,  the  digital  tapes  in  file  at  the  Lanont-Doher 
ty  Observatory  were  studied.  V.'e  have  selected  for  our  study  3 earth- 
quakes and  5 nuclear  explosions.  Of  those  only  37  records  were  found 
of  any  use.  Some  of  the  problems  found  were  that  the  records  were 
not  in  the  Lamont  file,  were  copies  defective,  failures  were  present 
in  the  digital  system,  and  magnification  was  too  low.  The  process 
of  finding  the  proper  place  in  the  topes  is  also  very  time  consu- 
ming,  since  there  is  not  a fast  search  system.  Outputs  of  the  digi- 
tal data  in  magnetic  tape  and  printed  form  ore  now  being  analyzed. 
Another  problem  encountered  is  using  HGLP-DIG'  data  is  the  variation 
of  the  ground  motion.  In  the  ALQ  (Albuquerque  Observatory)  station 
the  values  found  in  10  calibrations  from  Oct  73  to  Mar  74  gives  for 
the  mollification  at  30  sec  in  counts  per  micron  values  ranging 
from  302  to  1232  in  the  Z component,  5^2  to  852  in  NS  and  593  to 
965  in  EW  (A.  Murphy,  personal  comunication) . It  is  not  knwon  whe- 
ther the  changes  we re  tha  result  of  modifications  mode  to  the  system 
or  spontaneous  changes.  If  the  latter  is  true,  large  errors  con  be 
introduced  when  using  digital  data  from  stations  where  calibration 
is  not  repeated  so  often.  This  lock  of  stability  in  the  digital 
system  should  be  reflected  in  the  comparison  of  the  source  porome- 
ters  which  require  ground  motion  reduction  os  determined  from  ana- 
log and  digital  records. 


Determination  of  source  parameters 
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For  many  years  all  that  could  bo  found  about  the  mecha- 
nism of  earthquakes  was  the  orientation  of  two  -nodal  pianos  of  P 
waves,  one  of  them  corresponding  to  the  plane  of  fracture  'uid  the 
energy  released  as  seismic  waves  derived  usually  from  the  vaLe  of 
magnitude.  Progress  in  the  study  of  focal  mechanism  has  led  to  the 
defin'tion  and  determination  of  several  parameters  describing  more 
fully  the  processes  at  the  source  of  earthquakes  (Aki,  1972).  Those 
are  the  dimensions  of  the  pLnne  of  fracture,  the  seismic  moment, 
apparent  average  stress  and  stress  drop.  Of  the  several  methods 
proposed  for  the  determinations  of  these  parameters  we  hove  selec- 
ted two  shich  ore  mutually  independent,  one  based  on  the  spectral 
amplitudes  of  the  Rayleigh  waves  and  the  other  on  the  amplitude  and 
form  of  the  shear  wave  spectrum. 


Method  of  snectrol  amplitudes  of  Rsvlcinh  waves 

i i i»  ■ ■■.■■■■■■■  — ii  < 1 ■ ■ ■■  ■ ■ m n ii  * — — —■■■  i n —i  — 

Following  the  work  of  Ben-i-Ienahen,  .Jarosch  and  Roseman 
(1970)  the  Seismic  moment  con  be  expresed  in  terms  of  the  spectral 
amplitudes  of  the  vertical  component  of  the  Rayleigh  wave  in  the 
form 


Mo 


M UzM 


V^-t- 


(fRR+  sRSR«-i  ) 


(i) 


The  parameters  of  this  equation  ore  given  in  a convenient  fom  in 
tha  tables  of  the  same  article. 


The  length  of  the  fault  can  be  obtained  assuming  a propa- 
gating rupture  of  rectangular  for;,  from  the  "directivity  function". 


This  function  has  been  defined  by  Bcn-Mennhen  (1961)  as  the  ratio 
of  the  spectral  amplitudes  of  surface  waves  corresponding  to  rays 
leaving  the  focus  in  opposite  direction.  The  limitation  of  using 
only  data  from  stations  1 PQ  degrees  anart  ip  azimuth  is  not  neces 
sary  in  many  coses  (Udias,  1971).  If  o(  is  the  angle  between  the 
azimuths  of  two  stations  the  lenght  of  the  fault  is  given  by 


t = 


frnoot  (~  - CoS  (e+oP)  ) fm'm  (£*  ~ Cos  © ) 


(2) 


where  ffflax  and  fm^n  are  tne  frequencies  of  the  first  maximum  or 
minimum  of  the  directivity  function,  v the  ruoture  velocity,  CD  the 


phase  velocity  and 


the  azimuth  from  the  fault  clone . 


Method  of  shear-wave  soectrnm,  corner-frecmenc”' 

The  amplitude  end  form  of  the  spectrum  of  shear  waves  have 
been  related  to  the  values  of  certain  parameters  at  tne  source  (A hi. 
1972).  The  spectral  corner- frequency  fQ,  or  frequency  where  the  cue 
trum  bends  from  a constant  value  Oo  at  low  frequencies  to  a dimi- 
nishing value  with  frequency.  (Brune,  1970),  and  the  value  of  fl0 
itself  are  the  two  parameters  of  the  spectrum  related  with  source 
cncractoristi.cs,  A third  parameter  obtainable  from  the  spectrum  is 
€ wich  controls  the  high  frequency  decay  of  the  spectrum.  Accen- 
ting the  general  validity  of  Bruno's  model  and  the  physical  inter- 
pretation .jiven  by  him  to  the  three  parameters  fQ,  and  £ they 
can  be  used  in  the  specification  of  the  characteristics  of  the 
source  of  earthquakes  (Honk  and  Thatcher,  1 072). 
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The  spectral  amplitude  is  related  to  the  seismic 
momont  by  the  formula 

M0  = p (-7) 

The  corner  frequency  f is  related  to  the  radius  of  the  equivalent 
circular  fault  area  r by 

r - 2.34  (It ) 

2TI  f0 

The  stress  droo,  difference  in  shear  stress  acting  before  and  after 
an  earthquake  is  given  by 

A^-  ZMe  = i06pR^-°fo  (5) 


The  parameter  e is  related  to  the  ratio  of  the  stress  droo  ■AO’  to 
the  shear  stress  difference  on  the  fault  surface,  or  effective  she^r 
stress  Oipf 

c = ' (6) 

Oepf 


If  e is  taken  to  be  unity,  the  apparent  average  stress  defined  by 
Wyss  (1970)  as 
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rj  <r  = 21.1  p R-  Po  (8) 

This  formula  is  also  valid  to  a certain  approximation  for  the  cases 
where  6<1.  In  all  tnese  formulas  p>  is  the  shear  wave  velocity  a^d 
R the  reference  hypocentral  distance. 

Application  to  the  earthquake  of  6 June  1978 

Origin  time  and  locotion: 

6 June  1972,  05:25:50.2 

32.9  N , 39.9  "I,  h = 33 
°b  = 5.5 

Orientation  of  nodal  clones; 

plane  A:  N 66  W - 68  MV/ 
plane  B:  N 19  E - 76  NE 

Seismic  moment 

Seismic  moment  has  been  calculated  from  the  values  of  the 
smoothed  spectral  amplitudes  of  the  vertical  component  of  Eayleig' 
waves  corresponding  to  periods  of  50  and  100  seconds  and  from  the 
low-frequency  limit  of  the  amplitude  spectrum  of  shear  wav.  s , Re- 
sults are  listed  below: 
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Rayleigh  waves 


Station 

Mo 

dyne-cm 

T = 

50  sec 

T = 100 

ALQ  HGLP-ANL 

1025 

l\  102^ 

TLO  IIGLF-ANL 

3 

to21* 

3 I024 

KON  HGLP-ANL 

3 

1024 

.6  1024 

ALQ  HGLP -DIG 

3 

1024 

ALQ  Y/Y/5SN 

3 

1025 

102^ 

LPB  WV/SSN 

8 

102^ 

2 105 

Shear  waves 

i 

KON  HGLP-ANL 

5 1024 

ALQ  HGLP-DIG 

6.5  1024 

• 

Agreement  between  the  different  values  is  good.  As  an 

pi 

average  the  value  of  M is  3 lO^4  dyne-cra.  This  agrees  with  values 
obtained  for  the  moment  in  the  same  magnitude  range  by  other  au- 
thors (Aki,  1972).  The  agreement  between  the  values  obtained  form 
the  digital  and  analog  records  of  the  HGLP  is  also  good,  as  it  is 
betv/een  those  found  using  surface  and  shear  waves. 

In  figure  2 we  show  the  spectra  of  the  Rayleigh  v/aves 
recorded  at  the  same  station  (ALQ)  by  tie  V/Y/SSN , HGLP-ANL  and  1IGIP-DIC 
instruments,  reduced  to  grouse  motion.  The  spectrum  of  the  HGLP -PIG 
gives  smaller  amplitudes  for  perioas  shorter  than  26  sec. 

Length  of  foil  1 ting 

Length  of  faulting  has  been  calculated  from  the  directivity 
function  of  Rayleigh  waves  and  from  the  corner  frequency  of  shear 


r 


waves.  Results  aro  listed  below: 
Directivity 


Stations 

b km 

v km/sec 

KON/TLO  KGLP-ANL 

31 

1.5 

40 

2.1 

52 

3.0 

KON/ALQ  I1GLP-ANL 

46 

1 .0 

68 

1.5 

89 

2.0 

TLO/ALQ  HGLP-AKL 

42 

1.5 

52 

2.0 

70' 

3.0 

Corner  frequency 

r km 

KON  HGLP-ANL 

20 

ALQ  HGLP-DIG 

13 

The  values  from  the  directivity  depend  on  the  assumption 
of  a velocity  of  fracture  v.  In  previous  studios  (Udlas,  1 9^2 ) 
wo  have  seen  that  for  this  range  of  magnitudes  the  most  consistent 
value  of  v is  about  2 km/sec.  This  method  assumes  a rectangular 
fault  with  a propagating  fracture.  As  the  value  from  KNO/ALQ  is 
higher  than  the  other  two,  we  can  safely  take  b = 45  km.  From  the 
corner  frequency  v/e  obtain  the  radius  r of  the  equivalent  circular 
fault.  Taking  the  moan  of  the  two  values  r = 16.5;  this  an 
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an  arcs  of  855  km2  for  the  fault  surface.  This  area  in  the  rectan- 
gular model  of  lentgh  45  km  wil1  result  in  a fault  width  cf  19  km, 
which  is  also  a reasonable  value  for  a shallow  shock  of  this  mag- 
nitude. 

I 

Stress-drop  and  amvaront-stress 

From  the  corner  frequency  wc  can  also  obtain  values  for 
the  stress  drop  and  apparent  stress.  Using  the  HGLP  records  of  tv.'o 
stations  the  values  obtained  ore 

i 2 

KON  HGLP-AHL  Act  = 0.27  bars  ^cr  = 5 10  1 dyne/cm 

5 2 

ALQ  HGLP-DIG  AC  =4.1  bars  rfi-  3 10  dyne  /cm 

These  two  values  are  somewhat  low  specially  that  of  KON • The  value 
from  ALQ  io  nearer  to  what  is  usually  found.  The  aporent  stress 
can  be  found  directly  from  the  values  of  the  energy  and  the  moment, 
using  equation  (7). 

11  2 

Using  a value  for  y.  = 3.3  10  dyne/cm  and  for  the  ener- 
gy that  derived  from  the  magnitude  according  to  the  formula 

log  Es=  5.3  +,2.4  tn.Q 

which  gives 


the  resulting  value  for  the  apparent  stress  is 

= 1 . 1 10^  dyne/cm2 

This  value  agrees  better  with  other  values  of  this  parameter  for 
similar  magnitudes.  The  smaller  values  found  from  the  corner  fre- 
quency may  be  duo  in  part  to  the  fact  that  the  exprestens  used  as- 
sumed that  (■  - unity,  which  may  not  be  a realistic  assu  nation. 

Up  to  this  point  v;e  can  point  to  the  good  possibilities 
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offered  by  the  HGI.P  record  in  the  analysis  of  source  parameter  from 
shocks  in  the  magnitude  range  5.0^M<6;0.  So  fur  it  has  been  found 
also  that  outputs  of  the  digital  system  give  very  consistent  resits. 

P r opa  ga  t i on  v>a  ra  v e l e r s 

Two  other  areas  that  are  under  analysis  to  tost  the  reli- 
ability of  HGLP  data  are  the  determination  of  absortion  coeficie^t 
and  anelastic  rioperties  along  the  propagating  path  and  the  grout) 
and  phase  velocities  of  surface  waves.  Computer  programs  for  the 
processing  of  the  data  have  been  prepared.  These  involved  tie  use 
of  multiple  filter  techniques  to  determine  the  group  velocities  of 
the  energy  contained  in  narrow  bands  of  frequencies,  of  sliding- 
filters  to  eliminate  unwanted  disturbances  inthe  surface  wave  sig- 
nal and  programs  to  calculate  phase  velocities  by  different  methods. 

These  programs  will  be  applied  to  the  analysis  of  phase 
velocity  along  the  path  from  the  two  HGLP  stations  in  the  Medite- 
rranean TLO  and  EIL.  This  path  crosses  the  Mediterranean  sea  and 
it  is  exnedted  that  the  data  will  furnish  new  information  for  the 
interpretation  of  the  structure  of  thiscomplicnted  tectonical  area. 
Since  this  two  stations  are  also  W.V33N  stations  comparison  of  bob 
type  of  data  will  show  the  improve  .ent  of  the  results  due  to  the 
use  of  the  HGLP  stations. 
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